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Abstract. - The use of degree-degree correlations to model realistic networks which are charac- 
terized by their Pearson's coefficient, has become widespread. However the effect on how different 
correlation algorithms produce different results on processes on top of them, has not yet been 
discussed. In this letter, using different correlation algorithms to generate assortative networks, 
we show that for very assortative networks the behavior of the main observables in percolation 
processes depends on the algorithm used to build the network. The different alghoritms used here 
introduce different inner structures that are missed in Pearson's coefficient. We explain the differ- 
ent behaviors through a generalization of Pearson's coefficient that allows to study the correlations 
at chemical distances I from a root node. We apply our findings to real networks. 



Introduction. — In the last two decades, the use of 
complex networks in the study of many processes, such as 
the spread of diseases, random or intentional attacks, syn- 
chronization, etc. iJaiSj, has led to a deeper understanding 
\ of these processes. For example, the use of network-based 
models in epidemiology has demonstrated how the topol- 
, ogy affects the total fraction of the infected population 
' and how it's knowledge can be used to develop efhcient 
' immunization strategies [4]. Each of the aspects of the 
network topology has a substantial effect on the processes 
taking place on the underlying network. 
' A first feature of the network topology is the degree dis- 
tribution Pd{k), i.e., the fraction of nodes with degree k. 
One of the most used networks is the Erdos Renyi (ER) [S] 
with Pd{k) = e~^''^ (k)^ /k\, where (fc) is the average de- 
gree. However in real networks the degree distribution 
fits better to a scale-free (SF) Pd{k) ^ k^^ , where A is 
the broadness of the distribution. Many of these networks 
also have an exponential cutoff k [6J, and the degree dis- 
tribution can be represented by Pd{k) ^ k^^e^'^/'^ . 

Several analytical and numerical models on degree- 
degree uncorrelated networks were developed with the 
above distributions [ZllH]. However, in the last years it 
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has been observed that in real networks, the degrees of 
the linked nodes are correlated. For example, it is known 
that in social networks, nodes tend to be linked to others 
with similar degree, while in technological networks, large 
degree nodes tend to be connected with low-degree ones. 
Networks with the first pattern of connections are called 
assortative, while those with the second pattern are disas- 
sortative. Recent researches suggest that networks adopt 
a correlated structure in order to optimize some processes 
that are developed on top of them, such as synchroniza- 
tion 0, transport [TU], traffic dynamic [TT], congestion [T^ 
and growth [13]. As a consequence, it is necessary to con- 
sider models of networks with degree-degree correlations 
in order to study how the processes are affected by them. 
There are different measures to quantify the correlations 
between nearest neighbors [HIS]- The most detailed mea- 
sure is the joint degree distribution Pe{k,k'), which indi- 
cates the probability that a node of degree k is linked with 
another of degree k' . However, due to the large amount 
of information contained in this distribution, it is diffi- 
cult to interpret the results based on Pe{k,k'). Hence 
more global measures are needed to simplify the inter- 
pretation of the results. For example, a measure which 
follows from the above and is more simplified is the aver- 
age nearest neighbor degree of a node with connectivity k, 
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knn{k) — J2k' k'Pe{k\k'). The sign of the slope of knn{k) 
is positive (negative) for assortative (disassortative) net- 
works. However, the most used measurement is Pearson's 
coefficient, which is less detailed than the previous ones 
and gives a general overview of the correlation to nearest 
neighbors on the network. Pearson's coefficient is given by 

(fcfc%-((fc + fc')/2)^ 

((F+fc'2)/2)e-((fc + fc')/2)r ^' 

where e is the average over Pe{k, k'). Assortative networks 
have r > 0, while disassortative networks have r < 0. 
Finally, networks with r = are called uncorrelated. 

Many studies that concentrate on how the degree-degree 
correlation affects the processes on top of networks, use 
Pearson's coefficient to quantify the results. For example 
in Ref. [5] it has been observed that synchronization is 
enhanced for r ^ 0. Also, in Ref. [T^], under a gradient 
network formalism, the authors explain why in disassorta- 
tive networks the congestion is lower than in uncorrelated 
and assortative ones. On the other hand, some studies 
on the effect of degree-degree correlations were focused on 
the link percolation process [T4Ul7| . The investigations 
on this process are very useful since they provide infor- 
mation about how resistant a network is against random 
failures and also due to the relation between percolation 
with some aspects of the spread of diseases [T8H22]. For 
uncorrelated networks, it is well known that in percolation 
processes there is a critical fraction p = oi nodes/links, 
above which the network undergoes a second order phase 
transition f7[\W\ . The order parameter of the transition is 
Poo(p), which is the fraction of nodes that belong to the 
largest component for a given value of p. For infinite sys- 
tems (TV — >■ oo) and close to the threshold (pc), the order 
parameter behaves as [23] 

Pooip) -{p- Pcf, P>Pc + Sp, Sp > 0; (2) 

and the average value of the size of the finite clusters {s{p)) 
goes as 

(s(p)) - (p-pc)"'^, p<Pc-Sp, Sp> 0; (3) 

where P — I and 7 = f in the mean field approach |23j . 
In the case of finite networks with N nodes, Eq. (jS]) obey 
the general scaling law 

(s) ^ {sip,N))=N^^f[{p^p,{N))N''), (4) 

where / is a scaling function that behaves as a constant at 
the threshold, with (s) ~ 7V^/^ for uncorrelated ER and 
SF networks with A > 4 [241125]. 

Newman [26] and Miller [27] showed that networks with 
r > are more resilient to damage than those with r < 0. 
However, there has been little discussion in the literature 
about whether the percolation phase transition remains 
of second order depending on the value of r. This issue 
was only addressed in Ref. ^28j where they found that in 



correlated networks constructed using a Metropolis algo- 
rithm and characterized by r, for the disassortative case 
percolation belong to the same universality class as un- 
correlated networks but highly assortative networks have 
no second order phase transition. A recent investigation 
by Dorotgsev et. al [20] pointed out that recursive trees 
can be constructed with a highly correlated structure but 
with zero Pearson's coefficient. On the other hand, in 
Ref. [30] the authors show that positive correlations af- 
fects the loop statistics on networks. Thus, Pearson's co- 
efficient could hide a more complex pattern of connections 
between nodes that also affect the processes that spread 
on networks. This raises the question as to whether or 
not Pearson's coefficient is a reliable tool for characteriz- 
ing correlated networks and, in particular, the percolation 
transition on assortative networks. 

In this letter we use different correlation algorithms to 
determine the effect of very assortative networks on the 
link percolation transition. We find that the results de- 
pend strongly on the algorithm used to correlate the net- 
works and not on the value of Pearson's coefficient. More- 
over, we find that strong correlations above the first shell 
on assortative networks affect the percolation transition. 

Correlation Algorithms. — We will focus only on 
assortative networks because we want to understand the 
effect of different algorithms in the percolation transition. 

In this work we use rewiring algorithms to correlate 
the networks, since they preserve the degree distribution 
Pd{k). In particular we focus on two correlation algo- 
rithms: the Exponential Random Graph (ERG) |28| and 
the Local Optimal Algorithm (LOA) introduced by us. 

For the ERG algorithm, the process to correlate a net- 
work or graph G [2H], is based on a Metropolis dynamic, 
which uses a "Hamiltonian" function given by H(G) = 
—J Aijkikj, where Aij is the component of the ad- 

jacency matrix, ki (kj) is the connectivity of node i (j) 
and J is a control parameter for the correlation. For 
J > (J < 0) an assortative (disassortative) network 
is built. This algorithm generates a Gibbsian ensemble 
network and is ergodic. The algorithm correlates the net- 
work, successively applying the following steps: 

1. in a network G, two links are chosen randomly to be 
correlated. 

2. rewiring is allowed (from configuration G to G') with 
probability w = min{l, exp[-(iJ(G") - H{G))]}, dis- 
allowing autoloops and multiple connections. 

The rewiring process continues until the steady state, 
where for large enough system sizes, r depends only on 
J (not shown here). 

In the LOA algorithm, a network is correlated as follows: 

1. two links are randomly chosen from the network to 
be correlated. 

2. for assortative (disassortative) networks, we choose 
the configuration that increases (decreases) most 



p-2 



Effect of degree correlations above tlie first shell on the percolation transition 



Pearson's coefficient r without generating autoloops 
or multiple connections. 

This algorithm is not ergodic but achieves the desired cor- 
relation faster, since in contrast to the previous one, we 
do not have to wait for the equilibrium value in order to 
reach the desired value of r. 

Results and Discussion. 

Link percolation on assortative networks. In this sec- 
tion we study the effects of the correlation algorithms on 
the percolation transition in assortative networks. 
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Fig. 1: Plot of Poo(p) for ER networks with (fc) = 2, 
AT = 8 X 10'' and r = 0.78, correlated with LOA (sohd 
line) and ERG (dashed line, J = 1). In the inset we show 
an enlargement of the main plot close to the threshold. 
The simulations were done over 10'^ realizations. 

In FiglU we show Poo(p) for an ER network with fixed 
N and r for the ERG and LOA models. Even though the 
curves do not differ significantly between them, in the inset 
we can see that Poo(p), close to the threshold computed 
with different algorithms converge to zero with different 
slopes, suggesting that the transitions are different even 
if the same value of r is used in both algorithms. The 
behavior of Poo{p) for networks generated with the ERG 
algorithm seems to be nonsingular. 

It is well known that in the thermodynamic limit a sec- 
ond order phase transition, has a singularity in the deriva- 
tive of the order parameter at the threshold Pc- In Fig [2] 
we plot the derivative of Poo(p) with respect to p for dif- 
ferent networks sizes. From the plots, we can see that in 
strong positive correlated networks with the LOA algo- 
rithm, the derivative of Poo(p) exhibit a singularity as we 
increase the system size, characteristic of a second order 
phase transition. However for the ERG model the slope is 
smooth and does not depend on N . Therefore, we confirm 
that the transitions for the ERG and the LOA are not the 
same |28| . 

In order to verify the existence of a second order phase 
transition, we compute (s) for ER and different network 
sizes as shown in Figl31 We can see that with the LOA 
correlation model as N increases, the peak of (s) increases 
around the critical threshold Pc{N) as in a second-order 




Fig. 2: Plot of the derivative of Poo as a function of p 
for ER with (fc) = 2, r = 0.78, correlated with LOA (a) 
and ERG (b) for: TV = 8 x 10* (sohd hue), N ^ 16 x 
10"* (dashed line), TV = 32 x 10^ (dash-dotted line) and 
= 64 X 10^ (dotted line). In the insets we show an 
enlargement of the main plots close to the threshold. The 
simulations were done over lO'^ realizations. 



phase transition [23]. At p = pc (see Eq. g])) (s) N°-^^, 
suggesting that percolation in assortative networks gener- 
ated by the LOA belong to a different universality class 
from that in uncorrelated networks. As noted by Noh us- 
ing the ERG model [28j, (s) has no peak independent of 
the value of N, which still indicates that the transition is 
not of second order. 

We run all our simulations for SF networks with expo- 
nential cutoff and find similar behavior that for ER net- 
works. In Fig. m we show that the derivative of Poo(p) 
and (s) do not exhibit any singularity as we increase the 
network size for the ERG model, in contrast to the LOA 
model. 

Despite the discrepancies between the assortative net- 
works, we find that disassortative networks generated by 
both algorithms undergo a second-order phase transition 
with (s) ^ iV^/'^, which not only confirms that disassor- 
tative networks belong to the same universality class as 
uncorrelated ones, but also that the critical exponents are 
independent of the algorithm used to generate negative 
correlations. Consequently, we conclude that while Pear- 
son's coefficient is a good measure to uniquely characterize 
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Fig. 3: Plot of (s) as a function o p for the same networks 
as in FiglH LOA (a) and ERG (b). In the inset of (a) we 
plot the height of the peak of (s) as a function of N on 
a log-log scale. The dashed line represents a power law 
fitting from which we obtain 76 = 0.16 (see Eq. Q). 



the percolation transition on disassortatives and uncorre- 
lated networks, it fails to characterize the transition on 
very assortative networks. The discrepancies between the 
assortative networks generated with different algorithms 
must underlie on structural differences introduced by the 
algorithms, which cannot be explained by only Pearson's 
coefficient. Hence, we will show below the different effects 
that both algorithms introduce on the network topology. 

Correlations above the first shell. In order to explain 
the differences between assortative networks with the same 
value of r but correlated with different algorithms, we pro- 
pose a new measure called generalized Pearson's coefficient 
rg, which measures the correlation between the connectiv- 
ity of a root node and the neighbors at a chemical distance 
£, i.e. which are located in the i — th shell from the root . 
In this way, we denote the usual Pearson's coefficient r as 
ri, Pearson's coefficient to second neighbors as r2 and so 
forth. We define, 



re = 



{kk')t- {{k + k')/2)l 
{{k^ + k'^)/2)e-{{k + k')/2)j 



(5) 



Fig. 4: Derivative of Poo{p) for SF networks with A = 2, 
K = 10 and r = 0.55 correlated with LOA (a) and ERG 
(b) for: iV = 8 X lO'^ (solid line), = 16 x lO'' (dashed 
line), A^ = 32 X lO'' (dash-dotted line) and A^ = 64 x 10'' 
(dotted line). In the inset we plot (s) as a function of p 
for LOA (a) and ERG (b). 



a node of degree k' at the £ — i/i-shell. 




where (. . .)e is the average over Pi{k, k^), where Pe{k, k'^) 
is the probability that a node of degree k is connected to 



Fig. 5: Plot of ri as a function of the shell number I for 
LOA (□) and ERG (Q) in ER networks with (fc) = 2, 
r = ri = 0.78 and A^ = 10'*. The dashed lines are used 
only as a guide to the eye. In the inset, we show It as 
a function of the shell number I. As seen in the plot 
I2 = 1.21, /3 = 1.46, Ia ^ 1.81. Then the ERG model 
generate stronger long range correlations than the LOA 
model. The simulation were done over 10"* realizations. 
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In FigOwe show for ER networks as a function of ^ 
in assortative networks Q for ri > 0. We can see that ri is 
a decreasing function of i. Imposing high assortativity in 
the network, correlations above the first shell from a root 
node are generated. Thus, correlation will build groups of 
nodes connected with similar degrees at shells close to the 
root, forming groups of nodes of low degrees and others 
with high degree. The high degree groups will form a 
strong core very resilient to random failures. The stronger 
the correlations are above the first shell, the more resilient 
groups are, smoothing out the transition. 

In the inset of FiglS] we plot the factor li defined as 
the ratio between for ERG and for LOA networks. 
Then > 1 (/^ < 1) indicates that ERG networks have 
higher (lower) degree correlation at a distance £ than LOA 
networks. It is easy to see that for shells above £ — 1, 

increases with £ reaching a value of I4 = 1.8 for the 
fourth shelQ for N = 10''. Thus ERG networks gener- 
ate stronger correlation above shell £ = 1 than LOA net- 
works supporting our picture that high assortative ERG 
networks posses a stronger structure of groups, with nodes 
of similar connectivities, than in the LOA model. This 
introduces two aspects of the network that contributes to 
erase the second-order phase transition in the ERG model. 
First, ERG networks have more homogeneous groups of 
nodes with low connectivity that are easier to fragment, 
generating smaller clusters than in LOA assortative net- 
works. On the other hand, in ERG networks there are 
also more homogeneous groups of highly connected nodes 
that remain in the largest component through the perco- 
lation procesH, preventing the network from fragmenting 
into large finite clusters, and smoothing out the curve of 
Poo(p) near the threshold (see Fig 15]). Consequently, these 
two aspects induce a greater presence of small clusters, 
leading to a non-diverging (s). To determine if this be- 
havior is representative of real networks, we measure rg 
on two different positive correlated SF real graphs, the 
condensed matter (cond-mat, A « 1.60 and k fa 40) [55] 
and mathematics (A « 1.72 and k w 35) [34J coauthorship 
networks. In both cases, we obtain that positive corre- 
lations are non-negative only until the second shell with 
ri = r = 0.18, r2 = 0.01 for the cond-mat collaboration 
network and ri = r — 0.12, r2 — 0.03 for the mathematics 
collaboration network. Consequently, real networks posses 
mainly correlation to first neighbors (r) , in contrast to the 
stronger correlation structure above £ — 1 found in the 
theoretical networks, particularly for ERG, and therefore 
the real graphs analyzed here have a less defined struc- 
ture of groups of nodes connected with similar degrees, 
suggesting that there are constrains which prevent net- 



^The understanding of the disassortative behavior of goes be- 
yond the scope of this work. 

^For a SF network with A = 2, k = 10, r = 0.55 and N = 10"', 
we obtain that I2 = 2.11, however positive correlations only reach 
to the second shell. 

■'The same picture can be seen in onion networks |31I32| . However 
onion and assortativity are distinct properties as pointed out in |ijlj. 
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Fig. 6: Plot of the derivative of as a function of p close 
to the threshold for a) the cond-mat coauthorship network 
with N = 40421 and r = 0.18 and b) the mathematics 
coauthorship network with N — 391529 and r = 0.12. The 
solid lines correspond to percolation on the original data, 
and dashed and dotted lines correspond to LOA and ERG 
networks, respectively with the same degree distribution 
as the real data. We can see that all the curves exhibit 
a pronounced steep growth at the threshold characteristic 
of a second order phase transition. 



works from evolving towards extreme correlations. This 
can be seen, for example, on both coauthorship networks 
used here. In general, coauthorship networks are made 
of research groups composed by senior researchers with 
high connectivity and young researchers generally with 
low connectivity. On one hand, senior researchers from 
different research groups have some connections between 
themselves in order to increase the resources and impor- 
tance of their groups. Moreover, young researchers from 
the same group collaborate between themselves, because 
they generally work on similar projects. This explains the 
assortativity of the networks. On the other hand, in order 
to share their knowledges, senior researchers are highly 
connected to the young ones of their group, decreasing 
the assortativity. Thus the network evolves with an as- 
sortativity correlation, but with a not too high Pearson's 
coefficient in order to improve the functionality of the full 
coauthorship network. This is why real assortative collab- 
oration networks are not extreme assortativity correlated. 
Moreover, since this positive correlation decreases with 
£ as we can see in both real networks used here and in 
Router graph studied in Ref. |35) . the correlation above 
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the first shell is almost zero due to the small Pearson's 
coefficient, which leads to a critical behavior of percola- 
tion process on these networks. In Figl6] we plot for the 
real networks, the derivative of Poa{p) for the original data 
and after applying the LOA and the ERG models for net- 
works with the same degree distribution. As we can see, 
there is a steep growth of the derivative near criticality, 
which is equally sharp in the LOA and ERG networks, 
supporting our argument that in real networks the sec- 
ond order phase transition exists. In turn, we note that 
ERG and LOA networks with the same degree distribu- 
tion and Pearson's coefficient as real networks, have also 
a second order phase percolation transition. Our results 
suggests that Pearson's coefficient is a good indicator of 
the percolation transition order when the network is not 
strongly correlated above the first shell, but for high pos- 
itive degree-degree correlations, different algorithms may 
generate different behaviors in percolation and therefore, 
Pearson's coefficient is no longer useful to indicate the be- 
havior of the transition. As a consequence, it becomes 
crucial to use measures that take into account the inner 
structure of the network, such as rg, that as we show here 
is a better indicator of whether a network will undergo or 
not a critical phase percolation transition. 

In summary, we find that the second order percolation 
transition on theoretical networks depends on the inner 
structure imposed by the algorithms used to build them, 
and not only on Pearson's coefficient. This means that 
Pearson's coefficient hides a long-range correlation that 
could change dramatically the behavior on top of them. 
We propose a new magnitude r^, which allows to explain 
the discrepancies between percolation transition on assor- 
tative networks generated by different correlation algo- 
rithms. For the real networks analyzed in this letter, we 
find that they posses mainly first neighbors correlation, 
and thus the percolation transition is of second order. 

* * * 

This work was supported by UNMdP and FONCyT 
(PICT 0293/2008). The authors thanks Ana L. Pastore y 
Piontti and Joel C. Miller for useful discussions. 

REFERENCES 

[1] BoccALETTi S., Latora V., MoRENO Y., Chavez M. 

and Hwang D., Physics Reports , 424 (2006) 175. 
[2] DoROGOVTSEV S. N. and Mendes J. F. F., Evolution 

of Networks (Oxford University Press) 2003. 
[3] DoROGOVTSEV S. N., GoLTSEV A. V. and Mendes J. 

F. F., Rev. Mod. Phys. , 80 (2008) 1275. 
[4] Cohen R., Havlin S. and ben Avraham D., Phys. Rev. 

Lett. , 91 (2003) 247901. 
[5] Erdos p. and Renyi A., Publications Mathematicae , 6 

(1959) 290. 

[6] Amaral L. a. N., Scala A., Barthelemy M. and 
Stanley H. E., Proc. Natl. Acad. Set. USA , 97 (2000) 
11149. 



[7] Newman M. E. J., Strogatz S. H. and Watts D. J., 

Phys. Rev. E , 64 (2001) 026118. 
[8] Miller J. C, Phys. Rev. E , 76 (2007) 010101. 
[9] Sorrentino F., di Bernardo M., Cuellar G. H. and 

Boccaletti S., Physica D , 224 (2006) 123. 
[10] Xue Y.-H., Wang J., Li L., He D. and Hu B., Phys. 

Rev. E , 81 (2010) 037101. 
[11] Sun J.-T., Wang S.-J., Huang Z.-G. and Wang Y.- 

H., Physica A , 388 (2009) 3244. 
[12] Pastore y Piontti A. L., Braunstein L. A. and 

Maori P. A., Physics Letters A , 374 (2010) 4658. 
[13] Shao J., BuLDYREV S. V., Braunstein L. a., Havlin 

S. and Stanley H. E., Phys. Rev. E , 80 (2009) 036105. 
[14] Vazquez A. and Moreno Y., Phys. Rev. E , 67 (2003) 

015101. 

[15] Xulvi-Brunet R. and Sokolov I. M., Phys. Rev. E , 

70 (2004) 066102. 
[16] MoREiRA A. A., Andrade J. S., Herrmann H. J. and 

Indekeu J. O., Phys. Rev. Lett. , 102 (2009) 018701. 
]17] HooYBERGHS H., Van Schaeybroeck B., Moreira 

A. A., Andrade J. S., Herrmann H. J. and Indekeu 

J. O., Phys. Rev. E , 81 (2010) 011102. 
[18] Callaway D. S., Newman M. E. J., Strogatz S. H. 

and Watts D. J., Phys. Rev. Lett. , 85 (2000) 5468. 
[19] Cohen R., Erez K., ben Avraham D. and Havlin S., 

Phys. Rev. Lett. , 85 (2000) 4626. 
[20] Cohen R., ben Avraham D. and Havlin S., Physical 

Review E , 66 (2002) 036113. 
[21] Newman M. E. J., Physical Review E , 66 (2002) 016128. 
[22] Grassberger P., Math. Biosci. , 63 (1983) 157. 
[23] Stauffer D. and Aharony A., Introduction to percola- 
tion theory (Taylor & Francis) 1985. 
[24] BoRNHOLDT S. and Schuster H., (Editors) Handbook 

of graphs and networks (Wiley- VCH, Berlin) 2003. 
[25] Wu Z., Lagorio C, Braunstein L. A., Cohen R., 

Havlin S. and Stanley H. E., Physical Review E , 75 

(2007) 066110. 

[26] Newman M. E. J., Physical Review E , 67 (2003) 026126. 

[27] Miller J. C., Phys. Rev. E , 80 (2009) 020901. 

[28] NOH J., Physical Review E , 76 (2007) 02611. 

[29] DoROGOVTSEV S. N., Ferreira a. L., Goltsev a. V. 

and Mendes J. F. F., Phys. Rev. E , 81 (2010) 031135. 
[30] NoH J., The European Physical Journal B , 66 (2008) 

251. 

[31] Schneider C. M., Moreira A. A., Andrade J. S., 
Havlin S. and Herrmann H. J., Proc. Natl. Acad. Sci. 
USA , 3838 (2011) 108. 

[32] Wu Z.-X. and Holme P., Phys. Rev. E , 84 (2011) 
026106. 

[33] Newman M. E. J., Proc. Natl. Acad. Sci. USA , 98 
(2001) 404. 

[34] Palla G., Farkas I. J., Pollner P., Derenyi I. and 
VicsEK T., New J. Phys. , 10 (2008) 123026. 

[35] EcHENiQUE P., Gomez-Gardenes J., Moreno Y. and 
Vazquez A., Phys. Rev. E , 71 (2005) 035102. 



p-6 



